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We study the photoemission and optical conductiv- 
ity response of the strongly correlated metallic system 
Ca x Sri- x VOz. We find that the basic features of the transfer 
of spectral weight in photoemission experiments and the un- 
usual lineshape of the optical response can be understood by 
modeling the system with a one band Hubbard model close to 
the Mott-Hubbard transition. We present a detailed compari- 
son between the low frequency experimental data and the cor- 
responding theoretical predictions obtained within the LISA 
method that is exact in the limit of large lattice connectivity. 
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The understanding of the redistributions of the spec- 
tral weight that occurs in strongly correlated systems is 
a fascinating problem of condensed matter physics |[]-[| . 

Most transition-metal oxides have narrow conduction 
bands originated in partially filled 3c?— orbitals where cor- 
relation effects are significantly relevant. One of the clas- 
sical examples is the V2O3 compound and its derivatives 
which exhibit a Mott transition, that is, a metal-insulator 
transition (MIT) driven by electronic correlations in the 
absence of magnetic ordering . In a recent paper Fu- 
jimori et al. IJ examined the photoemission spectra of 
several transition metal oxydes demonstrating the dra- 
matic redistribution of spectral weight as consequence of 
electronic correlations in those compounds. Here, we con- 
sider another correlated system, Ca x Sri- x VOs, which 
displays a notable transfer of weight as a function of x in 
its photoemission spectra |?],|8| , and very unusual features 
in the distribution of its low frequency optical response. 
This system is particularly interesting as it allows for 
a systematic study of the evolution of the spectral line- 
shapes with x. We shall argue that the basic features 
observed in a recent photoemission study and in the 
new optical conductivity results that are reported here, 
can be explained in terms of the solution of a single band 
Hubbard model close to the MIT. We shall present a 
detailed comparison between the experimental data and 
the theoretical predictions of the model treated within 
the Local Impurity Selfconsistent Approximation (LISA) 
which is exact in the limit of large connectivity, or equiv- 
alently, large dimensionality J9|— |TT[] . 

The Ca x Sri_ x VC>3 compound has a perovskite-type 
structure where the nominal valency of vanadium ions is 
4+, which leaves only one 3d— electron per vanadium ion 



in the conduction band. Thus, in a first approximation, 
we shall assume that the system can be modeled by a 
half-filled single band Hubbard hamiltonian. The hope 
is that this simple model may capture the basic features 
of the low energy physics that is experimentally observed. 
The hamiltonian of the Hubbard model reads, 
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where summation over repeated spin indices is assumed. 

The LISA method jlO| is a dynamical mean field theory 
that is exact in the limit q — > 00, where q is the connectiv- 
ity of the lattice. Consequently, the hopping amplitude 
is rescaled as ta — > -4= to obtain a non-trivial model 
B. The lattice hamiltonian is then mapped onto an 
equivalent impurity problem which is supplemented by a 
selfconsistency condition. The detailed derivation of the 
resulting mean field equations is already given in great 
detail elsewhere ]L0|— p~2|| . We shall assume a semi-circular 
bare density of states p°(e) = (2/ttL») v /1 - {e/Df, with 
t = y, which is realized in a Bethe lattice and on a fully 
connected fully frustrated version of the model [p^|-p^|. 
This density of states has the same edge behavior as one 
obtained from a three dimensional tight binding model. 
As there is only one 3d— electron in the Ca x Sr\- x V03 
system, we set = that renders the model half-filled. 
We use the second order perturbative calculation (IPT) 
to iteratively solve the associated impurity problem jl3) . 
It has already been shown elsewhere that this method 
is in excellent agreement with the exact solutions that 
can be obtained using exact diagonalization or quantum 
Monte Carlo techniques 0. Still, we use the IPT calcu- 
lation since the exact methods either lack the resolution 
or cannot reach the low temperature limit that is relevant 
for the present discussion. 

The solution of the Hubbard model in the limit of large 
dimensions has been obtained recently . One of the 
basic results is the existence of a MIT at an intermediate 
value of the interaction U c m 3D (U c — 3.37 D whithin 
IPT) . As the transition is approached by increasing U/D, 
one finds that the density of states splits into a character- 
istic three-peak structure. There is a central quasiparti- 
cle peak at the Fermi energy of width e* F that narrows as 
the transition is approached indicating the enhancement 
of the effective mass. This feature is accompanied by the 
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incoherent lower and upper Hubbard bands that emerge 
at frequencies of the order of ±t£ and which acquire the 
spectral weight that is transferred from the lower fre- 
quencies as the central peak narrows Oiyj . 

From the theoretical standpoint, the Ca x Sri- x VOs 
compound is extremely interesting, since the chemical 
substitution of a Sr 2+ ion by Ca 2+ of the same valence 
modifies the V - O - V bond angle from - 180° to 160°, 
thus controlling the strength of the ratio U/D. Therefore 
for concentrations ranging from x = to 1 we go from a 
less correlated to a strongly correlated metal. 

We now turn to the comparison between the experi- 
mental data and the theoretical predictions of the model 
Hamiltonian. The strategy will be to extract the rele- 
vant model parameters, i.e. U and D, from the experi- 
mental data, and use them as sole input for the theoreti- 
cal model. We emphasize that, following this procedure, 
there are no free parameters left to adjust in the theoret- 
ical fit of the experimental data. 

We first consider the recent experimental photocmis- 
sion data of Inoue et al. @ which we reproduce in Figj^. 
The low energy photoemission spectra consist of two con- 
tributions: a broad feature at higher frequencies and a 
another one, much smaller and narrow, that occurs at 
the Fermi energy. As we increase the value of x, i.e., 
the correlation strength, we observe a systematic transfer 
of spectral weight between the two spectral features . 
The lineshape of this spectrum is very reminiscent of the 
results for the density of states that are obtained from 
the large dimensional approach to the Hubbard model 
close to the MIT point (see inset of Figjl]). Thus, we 
associate the higher frequency feature to the lower Hub- 
bard band (LHB) and the feature at the Fermi level to 
the low frequency quasiparticle contribution. We extract 
the model parameters noting that: i) the LHB features 
peak at approximately — ^; ii) the half- width of the LHB 
should approximately correspond to D; and in) the ratio 
of the quasiparticle peak to the LHB spectral weight is 
controlled by the ratio U/D in our model. 

Photoemission experiments probe the occupied part of 
the density of states. Accordingly, we have multiplied the 
theoretical spectrum by the Fermi distribution function 
of 80K, and then convoluted with a Gaussian function 
of width 0.3eV that corresponds to the finite experimen- 
tal resolution JtJ. In Fig|l] we display the results for the 
theoretical photoemission intensity. The estimated val- 
ues for U and D are shown in the figure. We have check 
that a good fit can be also obtained by small variation of 
the parameters. Although one expects that D would also 
systematically change with composition, we choose to set 
equal to constant leV. This choice not only provides a 
reasonable fit, but it further allows to clearly appreci- 
ate the systematic evolution of the model lineshape as a 
function of the interaction. 

It is apparent that the theory can successfully capture 
the basic features of the experimental data. We observe 



the transfer of spectral weight between the LHB and the 
quasiparticle peak as a function of the interaction. In 
particular, we note that the theory captures correctly the 
intermediate frequency range around which the transfer 
of weight occurs (indicated with a shaded region in the 
figure). These results are more remarkable in regard to 
the complete failure of the LDA method to capture even 
the most basic features of the data (inset Fig.|l|) M. 
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FIG. 1. Above: Experimental photoemission spectrum of 
Ca x Sri- x VO'4 at 80K (from Q with background from oxy- 
gen 2p band subtracted). The inset contains the theoretical 
density of states p(ui) obtained from IPT at T = for a value 
of the interaction ^ = 3 close to the MIT. Below: Theoret- 
ical density of states p(u>) for the model parameter U indi- 
cated in the figure and D = leV obtained by IPT. The data 
corresponds to the ui < part of p(u) convoluted with the 
experimental resolution. The inset shows the corresponding 
LDA results (from 

Although we find that the basic features of the transfer 
of spectral weight are very well captured by our model, 
upon a more detailed comparison one notes that some 
discrepancies remain. Especially in regard of the shape 
of the small quasiparticle peak and the lower energy edge 
of the LHB. The discrepancies in the shape of the quasi- 
particle contribution, in particular, the renormalization 
of the value of p(ui) at the Fermi level, as was pointed 
out in Jt],||], may be due to a basic feature of the large-d 
approach, namely, the locality of the self-energy. This is 
probably connected to the overestimation by a factor of 
~ 4 of the enhancement of the calculated mass renormal- 
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ization compared with that deduced from the measured 
electronic specific heat coefficient (^ exp ~ 9.2mJ/K 2 mol 
for CaVOs 17|, this value is about 3 times larger than 
the LDA prediction). This issue may possibly be resolved 
by a better treatment of the spatial fluctuations. Never- 
theless, we emphasize that the goal of the present paper 
is not the description of the lineshape in its full detail, 
but rather, the understanding of some basic features in 
the distribution and transfer of spectral weight by the 
effect of correlations. 

We now turn to the comparison of our model results for 
the optical conductivity jl8| with the new experimental 
data that we report here. The optical reflectivity mea- 
surement was carried out for CaVO^ which is the most 
strongly correlated case. A single crystal sample was 
grown by the floating zone method and specularly pol- 
ished. It was annealed in air at ~ A70K for 24 hours in 
order to make the oxygen concentration become stoichio- 
metric jl9| . The reflectivity spectrum was measured at 
room temperature in energy range from 0.03eV to 6.5eV, 
using a Fourier transform spectrometer (0.03eV — 0.5eV) 
and a grating monochromator (0.5eV — 6.5eV). To ob- 
tain the optical conductivity by Kramers-Kronig analy- 
sis, we supplemented our data with that of Kasuya et 
al. pc[ | obtained in a higher energy region. Also a high 
frequency asymptotic function of the form u>~ was used, 
and we extrapolated our data to the lowest frequencies 
using the Hagen- Rubens formula. The result is displayed 
in Fig.|[ It basically presents three contributions: a low 
frequency narrow peak which is the Drude part of the 
optical conductivity, a second small feature that appears 
at energies around 1.75eV and a third large contribution 
that peaks at S.beV. The Drude part at low frequencies 
corresponds to the metallic character of the compound. 
However, unlike in normal metals, the Drude contribu- 
tion is only a fraction of the total optical spectral weight. 
This is a realization of the original arguments of Kohn 
for the optical response of systems close to a MIT |2l| . 

We can readily, and rather intuitively, interpret the 
higher frequencies peak structure in terms of excitations 
between the three features present in the density of states 
(c/. upper inset of Fig.Q). The large contribution in the 
conductivity spectrum at oj ~ 3.5e]/ corresponds to the 
direct excitations from the lower to the upper Hubbard 
bands, therefore this contribution should peak at an en- 
ergy of the order of the Hubbard interaction U . Accord- 
ingly, the position of the peak very approximately corre- 
sponds to twice the energy of the maximum of the lower 
Hubbard band in the photoemission spectrum of Fig^. 
Also, we note that the large spectral weight of this feature 
corresponds to the fact that most of the single particle 
spectral weight is concentrated in the incoherent Hub- 
bard bands. On the other hand, we should also expect op- 
tical excitations from the LHB to the quasiparticle peak 
and from this peak to the upper Hubbard band. As the 
energy for these two types of excitations is approximately 



^, they should also contribute to the optical spectrum 
with a feature at frequencies ui ~ ¥ « 1.75eV from the 
estimation made for this parameter in the photoemission 
discussion. This is in fact observed in the experimental 
results. We can further understand the small spectral 
weight of this feature from the fact that the quasiparti- 
cle peak in the p(ui) has only a small fraction of the total 
spectral intensity. 
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FIG. 2. The experimental optical conductivity (t(uj) in 
units of [ficrn]" 1 for CaV0 3 at T = 300K (thin line). The 
mid-infrared spectra (bold line) obtained after the subtrac- 
tion of the background contribution (~ (to — 2.5eV) 3 ^ 2 ) and a 
lorentzian fit of Drude part at lower frequencies (thin dashed 
lines). The model optical conductivity (bold dashed line) for 

= 3.2 and D — leV obtained by IPT (arbitrary units). 

In Fig.^ we also present the optical spectrum obtained 
from the LISA method using the model parameters de- 
termined before (x — 1). We find that all the basic fea- 
tures of our previous qualitative discussion are actually 
realized in the solution of the model. For a clearer com- 
parison of the model prediction to the experimental data, 
we have subtracted in the latter the high frequency back- 
ground due to excitations from the oxygen band to the 
UHB. From a fit to the lower energy tail of the oxy- 
gen band we estimate this background contribution to 
start at to ~ 2.5eV. Since our model does not contain 
disorder, the Drude part is a S— function at the origin. 
Thus, also for the sake of a clear comparison, we have 
subtracted the Drude part of the experimental spectra 
by fitting a lorentzian function at low frequencies. The 
remaining part of the experimental spectra then corre- 
sponds to the mid-infrared features at ^ and U that we 
discussed above. As the results show, these features are 
well captured by our model at the precise positions and 
with qualitatively correct relative weights. As in the case 
of the photoemission data, the remaining differences in 
the lineshapes may be possibly due to limitations of the 
IPT or the election of a semicircular p°. 

A final quantity that is useful in the interpretation of 
optical conductivity data is the ratio of the Drude part 
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of the optical response to the total spectral weight. In 
an uncorrelated system, this ratio is one. The effect of 
the correlations that drive the system close to the MIT 
is to split the contributions to cr(u>) into a 5-function 
like part at low frequency that corresponds to the coher- 
ent quasiparticle conduction and a regular mid-infrared 
contribution describing the incoherent motion. Within 
our present mean field theory, this ratio is obtained as 
R = Z^y, where Z is the quasiparticle residue and 
(K a ) is the expectation vaue of the kinetic energy in the 
free system (U = 0) Js 22 1. From the experiment we find 
that R eX p ~ 0.3, which is in good agreement with the 
prediction of the theory R = 0.25 ± 0.03 (the uncertainty 
comes from the approximate determination of the model 
parameters U and D and the rapid variation of R close 
to the MIT). We note that 1/R is usually called the op- 
tical mass renormalization and that in the d — ► oo limit 
— = 1JZ. Thus, in our mean field theory these two 
quantities are exactly related by the renormalization of 
the kinetic energy by the interaction. 

An important observation is that the consistent pic- 
ture that emerges from our simultaneous discussion of 
photoemission and ac-conductivity gives further support 
to the assumption that a single band model may be used 
for the study of the Ca x Sr\- x VO^ system. This obser- 
vation is rather surprising, as the typical magnitude of 
crystal fields ~ Q.leV is unlikely to explain how the 3- fold 
degeneracy of the ^g-bands is effectively lifted. A more 
plausible explanation may follow from the consideration 
of a more realistic degenerate model with inter-orbital ex- 
change and correlation interaction. The solution of such 
a model poses extra technical challenges and work along 
this line is currently under way. 

Finally, we would like to briefly mention that the 
temperature dependence of the dc-resistivity follows a 
p ~ AT 2 behavior which demonstrates that the system 
is, as our model solution, a Fermi liquid. 

The satisfactory agreement between the theory and 
the experimental data demonstrates that the Hubbard 
model treated within a mean field theory that becomes 
exact in the limit of large lattice coordination is indeed 
a good strating point for the study of electronic sys- 
tems where the competition of localization and itiner- 
ancy play a crucial role. The LISA solution of the single 
band model allows us to understand the basic features 
that are experimentally observed in the spectroscopy of 
Ca x Sri- x VOs. We find that the transfer of weight ob- 
served in photoemission data and the unusual distribu- 
tion of spectral weigh in the optical response are related 
to the proximity of the system to a MIT point character- 
ized by the emergence of the small energy scale [ p"5[ . 
Thus we obtain within a single model with essentially no 
adjustable parameters, a consistent interpretation of the 
unusual features observed in both experiments. 
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